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Abstract. Transformation optics is used to design a gateway that can block 
electromagnetic waves but allows the passage of other entities. Our conceptual 
device has the advantage that it can be realized with simple materials and structural 
parameters and can have a reasonably wide bandwidth. In particular, we show that 
our system can be implemented by using a magnetic photonic crystal structure that 
employs a square ray of ferrite rods, and as the field response of ferrites can be tuned by 
external magnetic fields, we end up with an electromagnetic gateway that can be open 
or shut using external fields. The functionality is also robust against the positional 
disorder of the rods that made up the photonic crystal. 



PACS numbers: 41.20.Jb, 42.25.Fx, 42.25.Gy 
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1. Introduction 



Transformation optics PQ 121 [3] has paved the way for the development of optical devices 
that can realize functionalities that were thought to be possible only in science fictions 
[1]-[T1]. One such conceptual device that has attracted great public interest is a gateway 
that can block electromagnetic waves but that allows the passage of other entities. This 
device can be viewed as an implementation of a "hidden portal" mentioned in fictions 
[I5l flG] . However, the feasibility of such devices is limited by the very complex material 
parameters and the narrow bandwidth. Here, we show that gateway-type devices can 
actually be realized with simple parameters and they can have wider band widths such 
that the concept is closer to reality [IZ|-[22] than previously thought. The structure 
can be implemented by using the magnetic photonic crystal structures that are field 
tunable, resulting in an invisible electromagnetic gateway that can be open or shut 
using magnetic fields. 



2. A new "superscatterer" 



We start from a simple implementation of transformation optics. In Fig. [T|a), an 
object (colored in blue) is placed to the left side of a double negative medium (DNM) 
with e = = —1. Let the object be of permittivity Eq and permeability fiQ. The 
detailed shape and the length scale of the structure are shown in Fig. [l|a). From the 
viewpoint of transformation optics, the whole structure is optically equivalent to another 
object described in Fig. [l]^b) for far field observers. The equivalent permittivity and 
permeability tensors are e = Eqc and /2 = fioc with a constant tensor c. We introduce 
some parameters for ease of reference following the coordinates Xi, X2, yi, y2 and the 
angle a in Fig. [Ij Let A = X2 — Xi be the waist of the object in the x-direction, 
r 
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where (x, y) are the position coordinates of the equivalent transformation medium in 
Fig. [l](b). The inequalities for the definition of the parameter, p, come from three 
different coordinate transformations and the three regions are distinguished by the two 
dashed lines in Fig. [l|^b). 

It can be proved that (see Appendix A): 
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The coordinate transformation is based on the recent work on transformation optics 
with complementary medium (or folded geometries) pSl [21]. The original object is 
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Figure 1. A simple electromagnetic structure comprising of two 
homogenous regions (left panel) and an equivalent structure that gives the 
same scattering (right panel) . (a) The schematic plot for the scatterer that appears 
electromagnetically larger than its physical size. The red region is the DNM, the blue 
region is the original isotropic material that is to the left of the DNM, and the green 
region represents the air background, xi, X2, yi, j/2 and a are the detailed coordinates 
of the structure. The dashed lines in the region (x > 0) give the virtual boundary for 
the equivalent transformation medium, (b) The equivalent transformation medium is 
indicated by the yellow region, whose parameters are implemented by equation 
The dashed lines distinguish the three different regions for the definition of p. 

mapped into the equivalent transformation medium. The DNM together with its mirror 
image (the boundary is marked by the dashed hne in Fig. [l|^a)) of air form a pair of 
complementary media, in the sense the phase accumulated in one segment is exactly 
cancelled by another. A more detailed discussion of the optical property of this specific 
geometry and the associated coordinate transformation can be found in Appendix A. 
The structure in Fig. [l|^a) can be treated as a form of a scattering amplifier [25] because 
its scattering cross section can be much larger than its geometric cross section. The 
amplification effect originates from the excitation of surface plasmons in the DNM and 
air interface. 

As a concrete example, we suppose that the original isotropic material (the blue 
region in Fig. [l|^a)) has a large value of permittivity, eq = —10000 and /io = 1, which 
may be treated approximately as a perfect electric conductor (PEC). Figure |2](a) shows 
the scattering field pattern of the structure in Fig. [l](a). The plane wave is incident from 
the top to the bottom and has transverse electric (TE) polarization, for which the E field 
is along ^-direction. In this paper, we only consider the TE modes for simplicity. Note 
that the same idea works for transverse magnetic (TM) modes as well. The frequency 
is 2GHz and xi = yi = 0.1 m, X2 = y2 = 0.2 m and a = n/A. The structure behaves 
like an equivalent material with parameters described by equation (|2| to the far-field 
observers. Due to the large mismatch of the impedance of the equivalent material with 
the air background, we expect that the equivalent object (the yellow domain in Fig. 
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Figure 2. The scattering pattern for the scattering object and the equivalent 
scatterer. (a) The scattering pattern of the structure depicted in Fig. [T]^a). (b) The 
scattering pattern of the equivalent material but with the PEC replacing the material 
implemented by equation (j2]). 



[Tj^b)) should scatter like a PEC, and for that reason we choose for comparison in Fig. 
|2]^b) the scattering field pattern of a perfect conductor filling up the entire domain of the 
virtual object in Fig. [l](b). The similar far-field scattered patterns between Fig. 
and Fig. [2](b) confirm the strong scattering effect. The permittivity and permeability 
of the DNM are actually chosen as —1 + 0.0001 x i in the simulations to avoid numerical 
divergence problems [21]. The PEC-like scatters are always chosen as the materials of 
permittivity at —10000 and the permeability is taken to be 1 in this paper for simplicity. 



3. An invisible electromagnetic gateway 

The amplified scattering effect can be utilized to make an invisible gateway (TB] . Suppose 
that a PEC wall separates the whole space into two regions, the upper domain and the 
lower domain. If there are channels (or gateways) opened in the PEC wall, people in 
the two different spaces can communicate with each other, both physically and through 
EM waves. However, if we replace the doors with the above-described configuration at 
a specific frequency, the communication for that frequency will be blocked because the 
systems behave like PECs. The most amazing fact is that the channel is in fact physically 
empty. There is nothing but air in the channel so objects can "walk through" but the 
channel is blocked as perceived by the eye because light at the designated frequency 
cannot penetrate. Figure |3] is a schematic plot of such a gateway based on the above 
idea. With the same scale as the one in Fig. |2](a), we demonstrate the properties of 
such a gateway. We suppose that there is a line source located at (0.05m, 0.4m) with a 
frequency of 2GHz in the upper domain. Figure |4](a) shows that the waves cannot pass 
through the gateway and are excluded from the lower domain. However, without the 
DNM, the waves can propagate into the lower region as shown in Fig. |4]^b). 

We now consider the band- width of such a device. To be more concrete, we consider 



A simple route to a tunable electromagnetic gateway 



5 




Figure 3. The computation domains for the electromagnetic gateway. The 

yellow regions are the perfectly matched layers (PMLs), the red region is the DNM, 
the green region is air, the blue regions are PECs. The point source is located at above 
the gateway and it is denoted by a brown circle. The region between the DNM and 
the right PEC is the so called "invisible gateway" . 



the following dispersions for the DNM: 



1 



/(/+^rr (3) 



where fp = 2.828 GHz, T = 0.1 MHz, F = 1.5, /o = 1 GHz and 7 = 0.075 MHz. 
When / = 2 GHz, the relative permittivity and permeability return to — 1 + 0.0001 x i 
used in the above simulations. When / = 1.6 GHz, the impedances do not match at the 
interface of the air and the DNM while the refractive index of the DNM is about —1.76. 
Figure |4]^c) shows the electric field pattern while Fig. ^d) shows the field without the 
DNM. The DNM can reduce the penetration of the waves from the upper region when 
compared with Fig. |4]^c) and Fig. |4|^d). However, as n < —1, there is still a "slit" 
between the virtual boundary and the right PEC, which allows the waves to propagate 
somewhat into the lower region. If n = — 1, the virtual boundary is simply the mirror of 
the interface of the DNM and PEC on the left with a; = as the mirror. The position 
of the virtual boundary can be obtained heuristically from the image-forming principle. 
When / = 2.4 GHz, the impedances are still mismatched while the refractive index of 
the DNM is about —0.56. As the absolute value of n decreases, the virtual boundary of 
the image expands and there will be no passage for the waves to penetrate because the 
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virtual boundary overlaps with the PEC on the right-band side. For example, we plot 
the electric field pattern in Fig. |4|^e) and the case without the DNM in Fig. |4]^f). We 
find that the DNM can eliminate the penetration of the waves from the upper region. 
That means that the negative band of the DNM beyond / = 2 GHz (or — 1 < n < 0) 
is the working frequency of the designed gateway. However, we also find that the wave 
blocking effect will become weaker with higher frequencies (in this case, about 2.5GHz) 
in our simulations. As negative index media are intrinsically dispersive, the bandwidth 
has to be finite as the parameters will deviate progressively from those required by 
transformation optics. But the present gateway is shown to be relatively robust and 
has a broad operation bandwidth of about 20%, which can be regarded as a broadband 
device. To have a broader bandwidth, we can simply reduce the distance of the DNM 
and the right PEC to enhance the overlapping while the functionality of the gateway 
is sacrificed. As there are extensive designs of the DNMs at various wavelengths (both 
theoretically and experimentally) [I7]-[22], it would be reasonably feasible for the present 
gateway to be realized in the future. 

If the losses of the metamaterial are large (i.e., the imaginary parts of the parameters 
of the DNM are larger than 0.1), the bandwidth of the gateway will be small and the 
functionality will be compromised. The absorption of the DNM is the key difficulty for 
both the gateway described here as for other devices such as the perfect lens. All our 
simulation results are calculated using the COMSOL Multiphysics finite element-based 
electromagnetics solver. 



4. The implementation of a tunable electromagnetic gateway 

Very recently, it has been demonstrated that the DNM can be realized with a simple 
array of ferrite rods without any metallic components [26], which can be used to 
implement the present gateway. One of the unique merits of this magnetic metamaterial 
is that its optical properties are magnetically tunable. As such, the gateway can be 
manipulated using external magnetic fields. The magnetic metamaterial is a periodic 
square array of subwavelength ferrite rods with radii of r = 3.5mm and a lattice constant 
of a = 10mm. The permittivity is taken to be e = 25, and the permeability of the 

/ fir -ifiK \ 

i/J-K Air 



ferrite rods has the form fi 



with f^r = l+ and 

VO l) 

= (i^g-2THa^f)2-4n2f2 ; wherc is the damping coefficient, uo = 'JiHq is the resonance 
frequency with 71 being the gyromagnetic ratio; Hq is the sum of the external magnetic 
field and shape anisotropy field along the ^-direction, ujm = TiM^ is the characteristic 
frequency with Mg the saturation magnetization and it is taken as Mg = 1750(7, typical 
for single-crystal yttrium-iron-garnet (YIG). As the absorption of single-crystal YIG is 
extremely low, we can set ai = 0. 

The numerical simulations on the gateway under different conditions are performed 
by using the multiple scattering method Figure [sl shows the electric field intensity in 
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Figure 4. The functionalities of the gateway (with or without the DNM) at 
different frequencies, (a) The electric field pattern for the gateway at 2GHz. (b) 
The electric field pattern for the gateway without the DNM at 2GHz. (c) Same as (a) 
but for 1.6GHz. (d) Same as (b) but for 1.6GHz. (e) Same as (a) but for 2AGHz. 
(f) Same as (b) but for 2AGHz. 
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Figure 5. The implementation of the invisible gateway using magnetic 
photonic crystals. Here, we show the electric field intensity pattern for 
diffterent geometries working at different values of the external magnetic 

field, (a) Hq = 500Oe, (b) Hq = 4750e, (c) without magnetic metamaterials, (d) 
with magnetic metamaterial replaced by PEC. 

the logarithmic scale. The bulk PEC is replaced by a discrete system of the periodically 
arranged PEC rods (marked with dark green solid circles in the figure), and the ferrite 
rods are denoted as hollow circles. Under Hq = 500Oe and the frequency / = 2.55 GHz, 
the effective refractive index of the magnetic metamaterial is n = — l(£e// = f^eff = 
In Fig. |5](a), it can be seen that the electric field is excluded from the air channel, so that 
the open channel appears to be blocked to the eye at this frequency. This demonstrates 
that the gateway can be implemented by a photonic crystal type structure in which each 
element (each ferrite rod in the periodic array) is identical. By changing the external 
field to Hq = 4750e, the refractive index of the magnetic metamaterial becomes n = 1( 
Eeff = fJ^eff = !)• The electric field intensity under such conditions is shown in Fig. 
ISFb), which shows that the channel is open for EM wave passage, or the passage appears 
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to be open to the eye. Our results demonstrate that the gateway can be implemented by 
a very simple configuration and the effect is tunable by an external field. As shown in 
Appendix B, the effect is fairly robust to the disorder of the position of ferrite cylinders. 
Figures [sj^c) and ^d) illustrate the fact that the channel is electromagnetically open if 
the ferrite rods are replaced by air and by PEC rods. In addition to the simple geometry 
and weak absorption, the magnetic metamaterial has a reasonably broad bandwidth as 
a DNM. 



5. Conclusion 



In conclusion, we showed that a robust and tunable electromagnetic gateway can 
be constructed using simple material parameters. The idea is conceived through 
consideration of transformation optics and can be realized using a photonic crystal 
type structure. 
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Appendix A 

In this appendix, we give the detailed coordinate transformation to produce the folded 
geometry used in the text. Let us first consider the following coordinate transformation. 



X = { -x\ - xf < x' < 0, 



(A.l) 



y' = y, z' 



where 



frf. yz-y rf. yz-y \ y < y < -y„ 



(xf, x|) = { (xi, X2), -yi<y <yi , (A.2) 

This coordinate transformation maps the virtual space (x-space) described in Fig. 
[TJ^b) into the physical space (x'-space) described in Fig. [l](a). If the virtual space is 
a vacuum, we can obtain the corresponding material parameters in the physical space. 
It contains two kinds of materials in the vacuum background, one a double negative 
medium (DNM) with e = jj, = —1 (shown in red in Fig. [l](a)), the other an anisotropic 
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Figure Al. The folded geometry used in the text. (a)A obliquely incident TE 
plane wave in free space. (b)The structure with the material described by equation 



(A.4| and the DNM are invisible. 



medium (shown in blue in Fig. [T](a)) whose parameters are derived in the following 
section. 

We can easily obtain the Jacobian transformation matrix, 
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(A.3) 



from equations (A.l) and (A. 2) for the mapping between the blue region in Fig. [T]^a) 
and the yellow region in Fig. lj(b), and we can also prove that = r and ^ = p in 
the text. From transformation optics, we can obtain the material parameters for the 
blue region. 



P 
r 





2 
r 
1 
r 










(A.4) 



Figure |Al[ a) shows the propagation of the incident transverse electric (TE) 
polarization plane wave from the left to the right but at an angle of 45 degrees from 
the x-direction in free space. The frequency is 2GHz. We consider the object described 



in equation (A.4) together with an DNM domian. The scales of the structure in our 
simulations are the same as the ones used in the text. Figure [Al| ^b) shows the scattering 
pattern of such a structure and demonstrates that the structure itself is invisible to the 
far-field observers. 

A modification of the above folded geometry can lead to a variety of conceptual 
devices. For example, if the virtual space is not a vacuum but an object of permittivity 
Eg and permeability /i<j, then the parameters in the blue region should be a permittivity 
tensor, EgEp, and a permeability tensor, figj^p- The object in the blue region together 
with the DNM domain will have a scattering cross section that can be much larger 
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Figure A2. The scattering pattern of an equivalent scatterer and the real 
scatterer. (a) The scattering pattern of an object of permittivity Ss = 2.25. (b) The 
corresponding scattering pattern of the scattering amphfier. 



than its geometric cross section. Figure |A2[ a) shows the scattering pattern of an 
object with permittivity Ss = 2.25 and permeabihty fig = 1. Figure |A2](b) shows 
the scattering pattern of the structure including the object with permittivity tensor 
EgEp and permeabihty tensor figp'p and the DNM. The two similar far-field patterns in 



Figs. A2 a) and |A2[ b) confirm the enhanced scattering effect. In the text, we apply 
this structure in an inverse manner. We suppose that the object in the blue region is of 
permittivity Eq and permeability fiQ. Then, the whole structure is equivalent to another 
object with its permittivity and permeability tensors, e = Eqc and /I = fioc (described 
in the text). We can easily obtain equation ^ by performing an inverse transformation. 

Appendix B 

In this appendix, we demonstrate the robustness for the implementation of the gateway 
using the magnetic photonic crystal structure. We introduce disorder by displacing the 
ferrite rods from their original positions with 

Ar = (a — 2r) X s X r„, (B.l) 

where a is the lattice constant, r is the radii of the ferrite rods, r„ is a uniform random 
variable on the interval (—0.5,0.5) (assures no overlap of cylinders), s is the strength 
of disorder. Figure |Bl[ a) corresponds to Fig. [5|a) but with a disorder of s = 0.5 
introduced. The gateway remains blocked to the eye for this strength of disorder. 
When Hq changes to 475 Oe, the gateway is open again as shown in Fig. |Bl[ b). Figure 
Bl c) corresponds to Fig. [sj^a) but with a strong disorder of s = 1.0 introduced. The 



functionality of the gateway is now slightly compromised but the blocking is at least 
partially effective at such a strength of disorder. With Hq = 475 Oe, the gateway is 
open again as shown in Fig. |Bl[ d). The results indicate that the functionality is robust 
to structural perturbations. 
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Figure Bl. The robust effect from the disorder. The electric field intensity 
pattern for the gateway with different levels of disorder working at different 
external magnetic fields, (a) s = 0.5, Hq = 500 Oe; (b) s — 0.5, Ho — 475 Oe; (c) 

s^l.O, Ho = 500 Oe; (d) s 1.0, iJo = 475 Oe. 
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